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1 Résumé 

> 

OO Dans de nombreux systèmes biologiques on trouve des structures arborescentes et bifurquantes comme les poumons, les 
04 arbres, les fougères, les coquilles des escargots, le système vasculaire, etc. Tous ces systèmes se distinguent par une in- 
variance d'échelle et se comportent donc quasiment comme des systèmes auto-similaires. Est-ce par soucis d'optimisation 
I qu'une telle symétrie a été retenue par la nature? Quels sont les comportements dynamiques et acoustiques de tels 
. systèmes auto-similaires ? Afin de donner une réponse il faut tout d'abord comprendre le rôle de l'auto-similarité dans 
^—i les comportements dynamiques comme la propagation d'ondes. Faute d'approches définitives qui décrivent la physique 
, de systèmes auto-similaires, on introduit ici le modèle simple d'une chaîne linéaire quasi-continue correspondant à une 
■ distribution auto-similaire de ressorts qui lient des masses ponctuelles entre-elles. Dans notre modèle la densité d'énergie 
L| ' élastique est une fonction exactement auto-similaire. L'équation du mouvement fournit une version auto-similaire de 
. ^ ' l'opérateur de Laplace dont le spectre de valeurs propres (la relation de dispersion) est une fonction de type Weierstrass- 
Mandelbrot. Cette fonction est exactement auto-similaire et est de plus une fonction fractale non-différentiable. 
?H En outre, il s'avère que dans le cas limite d'un milieu continu notre Laplacien adopte la forme d'intégrales frac- 

■ - - 1 tionnelles. Dans ce cas la densité de vibrations se révèle être une loi de puissance pour les basses fréquences avec un 
exposant strictement positif qui indique l'annulation de la densité de vibrations pour une fréquence nulle. Pour de plus 
amples détails nous renvoyons le lecteur à notre article récent [2] . 

Mots-clés: Auto-similarité, fonctions auto-similaires, transformations affines, fonction de Weierstrass-Mandelbrot, fonc- 
tions fractales, loi de puissance, intégrales fractionnelles. 
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2 Introduction 



Ce fut dans les années soixante-dix du dernier siècle que Benoît Mandelbrot déclencha une révolution sci- 
entifique en développant sa Géométrie Fractale [1]. Toutefois, les orgigincs de cette discipline remontent au 
moins jusqu'au 19ème siècle [3]. Depuis lors, les applications des structures fractales se développent de plus 
en plus, citons par exemple les antennes fractales [4, 5]. En outre, des réseaux de type Sierpinski révèlent des 
caractéristiques vibratoires intéresantes [6]. 

Tandis que des systèmes fractals et auto-similiares sont déjà un sujet d'études bien établi dans de nombreux 
domaines de la physique, ils attirent aussi de plus en plus d'intérêt en mécanique analytique et dans les 
sciences de l'ingénieur. Il faut se rendre compte qu'il y a encore une large méconnaissance de l'impact de 
r auto-similarité en tant que symétrie sur les comportements physiques. Cependant quelques pas initiaux ont 
d'ores et déjà été effectués [6, 7, 8, 9, 10, 11], alors que le développement d'une mécanique fractale se fait 
attendre. Faute d'une théorie définitive, il est donc désirable d'établir des modèles de sytèmes auto-similaires 
qui soient à la fois pertinents et autant que possible simples à analyser. C'est donc l'objectif de cet article 
que de développer un tel modèle. 

Nous explorons les comportements vibratoires d'une chaîne linéaire ayant des interactions inter-particulaires 
harmoniques réalisant une distribution spatiale auto-similaire. Quelques études ont été consacrées aux chaînes 
linéaires avec aussi bien des propriétés fractales [8, 14, 15, 16, 17] que des propriétés en dépendance exponen- 
tielle [12]. 

Au contraire de ces articles consacrés aux grilles discrètes, nous analysons ici les comportements oscil- 
latoires dans une chaîne linéaire qui est non seulement fractale mais aussi exactement auto-similaire par 
rapport à ses interactions inter-particulaires. Il semble que notre approche pluridisciplinaire soit utile pour la 
modélisation de certains problèmes en mécanique des fluides, notamment en milieux turbulents [18, 19]. 

La démonstration est développée comme suit : le § 3 est consacré à la construction de fonctions et 
opérateurs auto-similaires. Avec cette machinerie on déduit une sorte de Laplacien auto-similaire. Bien que 
la construction d'un tel opérateur ne soit pas unique, on démontre pourtant que cette manière de construction 
dispose d'une certaine justification physique. Dans une approximation continue notre Laplacien revêt la forme 
d'intégrales fractionnelles. 

Dans le § 4 on établit un modèle physique simple de chaîne linéaire de longueur infinie comportant 
des interactions harmoniques (ressorts) auto-similaires. L'auto-similarité requiert une distribution continue et 
homogène des particules. L'équation de mouvement apparaît sous forme d'une équation d'ondes auto-similaire 
qui contient notre Laplacien auto-similaire. Les valeurs propres de ce Laplacien fournissent la rélation de 
dispersion qui est une fonction de type Weierstrass-Mandelbrot qui comporte aussi bien l'auto-similarité que 
des propriétés fractales. 

3 Définition du problème afîin 

Tout d'abord nous définissons la notion d' "auto-similarité" en termes de fonctions et opérateurs. Nous 
dénommons une fonction scalaire ^(/i) exactemant auto-similaire par rapport à la variable h si la condition 

^{Nh) = k(t>{h) (1) 

est satisfaite pour toutes valeurs h > 0. Nous dénommons (1) le "problème afïin"^. est un paramètre 
fixe pour un problème donné et A = un ensemble continu de valeurs propres dont la bande admissible de 
l'exposant ô = est à déterminer. Une fonction qui satisfait (1) pour un certain N et une bande A = 
représente une solution (f)N,5{h) du problème (1) dont le paramètre fixe N est donné. Il va de soi que si une 
fonction çi)(/i) satisfait (1) elle satisfait aussi (f){N^h) = K^(j){h) oiî s € Z (zéro inclus). En choisissant s = —1, 
il est clair que nous récupérons le même problème en remplaçant A et A par et A~^ dans (1). Il s'ensuit 
que nous pouvons restreindre notre analyse aux valeurs N > 1 sans aucune perte de généralité. De plus, il 
est justifie pour des raisons physiques d'admettre A, A^ comme nombres réels et positifs (A, A G IR"*"). 

Pour aborder le problème affin (1) comme problème de valeurs propres, il convient d'introduire l'opérateur 
affin An comme suit 

^nous nous restreignons aux transformations affines de la forme h' = Nh + c with c = 0. 
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ANf{h) =: fiNh) (2) 
On vérifie facilement que l'opérateur affin est linéaire, c'est à dire que 

Ân (cihih) + C2h{h)) = ciMNh) + C2h{Nh) (3) 

et encore 

Â%f{h)= f{N'h), s = 0±l,±2,...±oo (4) 

Avec cela on peut définir des fonctions d'opérateurs. Soit gij) une fonction scalaire assez lisse pour être 
développée en série de McLaurin comme suit 

oo 

5(r) = 5:a,r^ (5) 
nous définissons une fonction de l'opérateur afïin correspondant à g par la série 

oo 

a{UN) = Y.''seÂ% (6) 

s=0 

oii ^ représente un paramètre scalaire. Cette fonction de l'operérateur afïin agit sur une fonction f{h) scalaire 
comme suit 

oo 

9{iÀN)f{h) = Y,cisi'f{N'h) (7) 

s=0 

La convergence de cette série est à vérifier pour toute fonction f{h) donnée. Il est aussi utile d'exprimer 
l'opérateur afïin en utilisation la notation explicite 



In TV 
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^^(/i) = e ^îkThT (8) 

Par ce formalisme nous sommes à même de construire les fonctions et opérateurs auto-similaires nécessaires 
à l'analyse, dans la § 4, du problème physique d'une chaîne linéaire auto-similaire. 

3.1 Construction de fonctions et opérateurs auto-similaires 

Une solution du problème afïin (1) s'écrit 

oo oo 

m = E ^-'Âhfih) = E ^-'fiN'h) (9) 

s=— oo s=— oo 

pour toute fonction f{h) dont la série (9) converge régulièrement pour tous h. On introduit l'opérateur 
auto-similaire comme suit 

oo 

= ^ A-^Â% (10) 

s=— oo 

qui satisfait la condition d'auto-similarité À^Tn = ATV. Nous soulignons le fait que nous pouvons nous 
restreindre aux cas > 1 (A, A G IR) et exclure le cas pathologique A^ = 1. 

Il est important de préciser pour quelles fonctions f{t) la série (9) est convergente. On trouve que pour 
qu'une fonction f{t) soit admissible, elle doit satisfaire les conditions {t > 0)[2] 



et 



\imf{t) = aoe (11) 



\im f{t) = Coot'^ (12) 

t— >oo 
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où ao,Coo sont des constantes. Les deux exposants o;,/3 G R peuvent revêtir des valeurs réelles, positives 
ou négatives avec (3 < a. Pour une telle fonction, la fonction de (9) est convergente et donc existe si 
l'exposant ô = est constraint par 

P<ô='^<a (13) 
Le cas /3 = comprend des fonctions bornées \f{t)\ < M ci-inclus quelques fonctions périodiques. 

3.2 Une variante auto-similaire de l'opérateur de Laplace 

Dans l'esprit de (9) et (10) on construit une fonction exactement auto-similaire telle que 

4>{x, h) = TNih) {u{x + h) + u{x -h)- 2u{x)) (14) 

où u{..) signifie un champ continu arbitraire et suffissamment lisse, et la dépendence de T]\f{h) en variable h 
exprime que l'opérateur afiin Aj^[h) n'agit que sur h, de manière que Aj^{h)v{x,h) = v{x,Nh). En posant 
^ = on a l'expression 

oo 

^{x,h)= Y. {u{x + N'h)+u{x- N'h)-2u{x)] (15) 

s=— oo 

laquelle est auto-similaire par rapport à la variable h, soit A]\f{h)(j){x, h) = (p{x,Nh) = ^~^(p{x,h). Par contre 
elle est une fonction régulière par rapport à la variable x. De plus, on exige que le champ u{x) soit Fourier 
transformable, ce qui impose comme condition que l'intégrale 

/oo 
\u{x) \ dx < OO (16) 
-oo 

existe. Il est utile de noter que 

u{x + h) + u{x -h)- 2u{x) = 4sinh2 f-^-) uix) = h'^-^u{x) + 0{ h'^) (17) 



2 dx J dx'^ 

en utilisant u{x ± h) = e^'^'^u{x) avec u{x + h) + u{x — h) — 2u{x) = (^e^d^ + e~^d^ — 2^ u{x). Ce qui est 

possible puisque u{x) est supposée suffissamment lisse. De (16) et (17) il suit que la série (15) converge si 
l'exposant 5 se retrouve dans l'intervalle 

0<*=-^<2 (18) 

Alors on introduit une variante du Laplacien auto-similaire de sorte que (avec i = N-^) [2] 

\5,N,h)<^) =■■ N'^H^, N-^h) = cf>ix, h) (19) 
où (f){x,h) est la série (15). En utilisant (17), le Laplacien (19) se réécrit comme une série d'opérateurs 

A(,.,.) = iMh) sinh^ (^^) = 4 f: iV-^^ sinh^ A) (20) 



Notamment on observe que la condition d'auto-similarité 

^{S,N,Nh) = ^^^{S,N,h) 

est en fait réspectée par notre Laplacien. 



^(S,N,Nh) - ^^^(S,N,h) (21) 
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3.3 Approximation continue - lien avec les intégrales fractionnelles 

Pour une évaluation numérique, il conviendrait d'employer une approximation continue pour le Laplacien 
(20) . Pour cela on pose N = 1 + e (avec < e <C 1 de sorte que e !^ In iV) oii e est un paramètre supposé 
petit. On définit de phis se = v avec du f« e et = (1 + e)~ ^ e^. Dans cette approche Ri é" devient une 
variable quasiment continue dont s parcourt tous les s G Z. (9) s'écrit donc 

"° 1 r°° 

(j)(h) = N-'^fiN'h) ^- e-^''f{hé')àv (22) 

en posant he^ = r (/i > 0) et ^ = dî; avec t{v — oo) = et t{v —>■ oo) = oo, on trouve que 

m-^-l ^dr (23) 
L'application de cette relation au Laplacien (20) donne 



/i^" P {ujx - t) + ujx + t) - 2uix)) 
(x)«-y^ r^S dr (24) 



Cette intégrale existe pour P < < ô < 2 avec /3 < — 1 du fait des conditions d'existence de l'intégrale 
(16) et de la relation (17). En effectuant deux intégrations partielles entre les limites r = et r = oo on 
réécrit (24), à condition que < ô < 2, sous la forme d'une convolution 

^9i\x-r\)-^iT)dT (25) 

avec le Laplacien unidimensionnel ^ et le noyau 

et avec g{\x\) = — Mn|a;| dans le cas (5 = 1. De plus, il conviendrait d'exprimer (25) en termes d'intégrales 
fractionnelles soit 

A(^=2-D,e,/^)^(X) « ^ -l)iD-2) (^"^'^ ^ (-l)""^-^.) ^M^) (27) 

OÙ Aiu{x) = -^u{x) dénomme encore le Laplacien unidimensionnel et on a introduit D = 2 — ô > qui 
est positif dans l'intervalle admissible < ô < 2. Nous verrons par la suite que D peut être définie dans 
l'intervalle < ô < 1 comme la dimension fractale estimée [20] de la courbe de la relation de dispersion 
correspondant à notre Laplacien. 

Dans (27) on avait introduit l'intégrale fractionnelle de Riemann-Liouville 'D~^ laquelle est définie par 
(voir par exemple [21, 22]) 



(28) 



oii T(D) indique la fonction F qui représente une généralisation de la fonction factorielle aux nombres non- 
entiers positifs D > 0. La fonction F est définie par 

T{D) = r T^-^e-^àT , D>0 (29) 
Jo 

Pour des entiers positifs D > 0, la, fonction de F reproduit la fonction factorielle F(Z)) = {D — 1)! avec 
D = l,2,..oo. 
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4 Le modèle physique 

On considère une chaîne linéaire d'une distribution quasi-continue de particules. Chaque point spatial x 
correspond à un "point matériel" c'est à dire une particule. La densité de masse supposée homogène est égale 
à un à chaque point de masse. Chaque particule est associée à un seul degré de liberté représenté par le champ 
u{x, t) où X est la coordonnée Lagragienne et t le temps. Chaque particule localisée à x est non-localcmcnt 
connectée par des ressorts d'intensité Ç'* (Ç = N~^, N G U > l,ô > 0) aux autres particules qui se trouvent 
aux positions x ± N^h {h > 0), avec s = 0, ±1, ±2, .. ib oo. Le Hamiltonien d'une telle chaîne s'écrit 

H = - (u^{x,t) + Vix,t,h)) dx (30) 
Dans l'esprit de (9) l'énergie élastique V{x, t, h) est construite de manière auto-similaire, soit ^ 

V(x,t,h) = ^fNih) {u{x,t)-u{x + h,t)f + (u{x,t)-u{x-h,t)f (31) 
oii TN{h) est l'opérateur auto-similaire de (10). On arrive donc avec ^ = = à 

oo 

V{x,t,h) = - i'\{u{x,t)-u{x + hN',t)f + {u{x,t)-u{x-hN',t)f'\ (32) 

s=— oo 

L'énergie élastique remplit donc la condition d'auto-similarité par rapport à /i, soit 

ijv(/i)V(x, t, h) = V{x, t, Nh) = C^V{x, t, h) (33) 

L'exigence que l'énergie élastique soit finie entraîne la nécessité de convergence de la série (32). Pour cela 
on trouve pour l'exposant 5 la contrainte [2] 

< (5 < 2 (34) 

Afin d'établir (34) on a utilisé également le fait que u{x,t) est Fourier-transformablc. L'inéquation (34) 
détermine l'intervalle oii l'énergie élastique (32) converge. L'équation de mouvement s'obtient par 



(oii ô./ôu indique la dérivée fonctionnelle) soit 



(92 7, °° 

= - Yl {Mx,t) -u{x + hN',t) -u{x-hN',t)) (36) 

s=— oo 

-Q^ = ^(ô,N,h)u{x,t} (37) 

avec le Laplacien auto-similaire A(j^^ /j-) de l'équation (20). Comme (32), l'équation du mouvement (36) existe 
si ô est dans l'intervalle (34). On peut réécrire (37) sous la forme compacte d'une équation d'ondes 

Oi5,N,h)u{x,t) = (38) 
oîi □(^^iVj/i) 6st V auto- similaire analogue de l'opérateur d'Alembertien d'ondes ayant la forme 

°(5,7V,h) = \8,N,h) - (39) 

Le d'Alembertien (39) avec le Laplacien (20) décrit la propagation d'ondes dans la chaîne auto-similaire dont 
le Hamiltonien est donné par (30) . Il semble que cette approche puisse constituer le point de départ pour une 
théorie générale de la propagation d'ondes dans des milieux auto-similaires. 



^Le facteur supplémentaire 1/2 évite un double comptage. 
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Le but suivant est de déterminer la relation de dispersion, laquelle se constitue par les valeurs pro- 
pres (négatives) du Laplacien semi-négatif défini (20). On utilise le fait que le champ u{x,t) est Fourier- 
transformable et que la fonction exponentielle e**^^ est fonction propre du Laplacien (20). En écrivant 

1 fOO 

u{x,t) = — ù{k,t)é^''dk (40) 

27r J-oo 

on se rend compte que les amplitudes de Fourier û{k,t) satisfont 

-^{k,t) = -Co\k)ù{k,t) (41) 

On obtient donc 

u^{kh)=4 f; iV-^^sin2(^^) (42) 

s=— oo 

La série (42) décrit une fonction de type Weierstrass-Mandelbrot qui est continue et qui est pour < 
S < l non-différentiable [1 , 20] . La fonction de Weierstrass-Mandelbrot satisfait toujours la condition d'auto- 
similarité (42), soit 

u?{Nkh) = N^uj'^{kh) (43) 

sur la totalité de son intervalle de convergence (34). 

On souligne le fait que seuls les exposants 6 dans l'intervalle (34) sont admissible pour l'Hamiltonien (30) 
pour définir un problème bien posé. Il fut démontré par Hardi [20] que pour (^N > 1 avec ^ = < 1, c'est 
à dire 

< (5 < 1 (44) 

la fonction de Weierstrass-Mandelbrot (42) est aussi une courbe fractale ayant une dimension fractale (di- 
mension de Hausdorfi") D = 2 — ô > 1. Les figures 2-4 montrent des courbes de dispersion pour des valeurs 
décroissantes de ô dans l'intervalle < 5 < 1. Il est bien visible que la décroissance de ô produit une croissance 
de la dimension fractale D. La hausse de la dimension fractale de la figure 2 à la figure 4 est bien visible par 
leur comportement de plus en plus turbulent. Dans la figure 4, la dimension fractale est D = 1.9 déjà près de 
la dimension 2 du plan. Par contre, la figure 1 correspond avec ô = 1.2 k un cas non-fractal. Afin d'évaluer 
(42) approximativement il convient de remplacer la série par une intégrale en utilisant la même approximation 
qu'au paragraphe 3.3 (e ~ InA^). Dans le cadre de cette approximation on trouve pour \k\h suffisamment 
"petit" {h > 0), c'est à dire dans le cas limite d'ondes longues 

u;\kh) ^ (45) 

oii (|/c|/i)^ doit être de l'ordre de grandeur de e ou plus petit. La relation de dispersion est donc caractérisée 
dans le cas limite d'ondes longues, par une loi de puissance de la forme de LD{k) w Const Ikf^"^. La constante 
C qui est introduite dans (45) est donnée par l'intégrale 

C = 2 (1-cosr) 



r°° 1 - cosr) , 



laquelle existe uniquement si ô est dans l'intervalle (34). Dans ce cas limite, on obtient la densité d'oscillateurs 
par [2] 

pH = 2ifl (47) 

ivr duj 

qui est normalisé de sorte que p{uj)duj compte le nombre d'oscillateurs, c'est à dire le nombre d'ondes qui se 
propagagent avec une certaine fréquence, laquelle se trouve dans l'intervalle [a;,a;-|-da;]. On obtient également 
pour la densité d'oscillateurs une loi de puissance sous la forme de 
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avec ô toujours dans l'intervalle (34). On observe donc que la puissance | — 1 est restreinte à l'intervalle 
< I — 1 < oo pour < S < 2. En particulier, on trouve l'annulation de la densité d'oscillateurs pour une 
fréquence nulle p{uj — > 0) = 0. 



5 Conclusions 

Nous avons construit des fonctions et opérateurs de manière auto-similaire en utilisant une certaine catégorie de 
fonctions traditionnelles admissibles. Avec cette approche il est possible de construire une sorte de Laplacien 
auto-similaire et par conséquent de l'opérateur d'Alembertien d'ondes. A l'aide de cette machinerie on déduit 
une équation de mouvement qui décrit la propagation d'ondes dans une chaîne auto-similaire. Nous croyons 
que ce-modèle est l'un des plus simple pour prendre en compte l'aspect d'auto-similarité d'un milieux. Nous 
démontrons également la correspondance entre notre approche et le calcul fractionnel qui fut employé plus 
tôt [9] afin d'aborder des problèmes physiques dans des milieux fractals. 

Nous espérons que notre approche pourra servir de point de départ pour développer une théorie de la 
propagation d'ondes dans des milieux auto-similaires et fractals dans des contextes pluridisciplinaires. 
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